Inversion of two-band superconductivity at the critical electron doping of (Mg,Al)B 2 
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Electron energy-loss spectroscopy (EELS) was combined with heat capacity measurements to 
follow the change of superconductivity with systematic Al doping of MgB2 . By using x-ray diffraction 
and Vegard's law to assess the actual Al content in the samples, changes in behavior were found to 
be much more in agreement with theoretical predictions than in earlier studies. EELS data show 
that cr-band hole states disappear above 33% Al, approximately the composition at which the a 
band Fermi surface is predicted to lose its cylindrical shape in reciprocal space and break apart into 
ellipsoidal pockets. At this composition, the a gap obtained from the heat capacity data falls to the 
level of the it gap, implying that band filling results in the loss of strong superconductivity on the 
a band. However, superconductivity is not quenched completely, but persists with T c < 7 K up to 
about 55% Al, the Al concentration at which the entire a band is predicted to fall below the Fermi 
surface. Since, in the region 0.33 < x < 0.55, only the ir band has appreciable density of states, it 
becomes the stronger of the 2 bands, thus inverting the 2-band hierarchy of MgB2. 



Although the two-band superconductivity was first 
studied theoretically in 1959 lj, its detailed experimen- 
tal study has not been possible until the discovery of 
superconductivity near 40 K in MgB 2 in 2001 2]. The 
high critical temperature T c results from two key ingredi- 
ents: First, the boron a band, which is deep beneath the 
Fermi energy Ep in graphite, sits just above the Fermi en- 
ergy and provides a high density of p xy hole states 0, Q . 
Second, the boron-boron Ei g bond oscillation is renor- 
malized downward in energy (softened), into the range 
where the phonon coupling to the hole states is nearly 
perfect 0, Q . The 2-dimensional a band carries along 
with it superconductivity in a 3-dimensional tt band with 
p z electron states, which has a smaller superconducting 
gap energy and weaker electron-phonon coupling [4|. 

Since the a band appears to be so crucial to supercon- 
ductivity, many experiments have attacked the central 
question of how superconductivity changes upon substi- 
tuting Al or C for Mg or B, respectively, thereby provid- 
ing extra electrons that fill its hole states. A common 
observation has been that such dopants destroy super- 
conductivity 0, S 0, E3, EH For Al doping, T c (x) falls 
with increasing Al concentration x, sometimes increas- 
ing slope at x » 0.3 and reaching zero for x — 0.5 to 
0.6 [a IS Eii ■ The loss of superconductivity at x f» 0.5 
correlates with changes seen in NMR data [l4_ [l5jl- op- 
tical reflectance 0], and infrared absorption @- Heat 
capacity measurements show an approximately lin- 
ear decrease of both gaps out to x = 0.4, suggesting that 
interband scattering does not affect these trends strongly. 
A rise of the Fermi level consistent with band filling has 
been inferred from x-ray absorption spectroscopy (XAS) 
16], but quantitative data is not available. 

Theoretical predictions are generally consistent with 
these observations, but also add unusual details. Elec- 
tronic structure calculations within the virtual crystal 
approximation predict that the cylindrical Fermi surfaces 
attributed to the boron a band collapse at the T point in 
reciprocal space when x ~ 0.3, separating the cylinders 
into strings of ellipsoidal pockets [lg ■ Complete filling of 



the o~ band holes was predicted to occur at x = 0.56. The 
topological change at x = 0.3 suggests abrupt changes in 
the properties occur for somewhat lower Al content than 
observed in experiments. One such change predicted to 
occur within the 2-band Eliashberg formalism 0, [2(j 
is a crossover in the hierarchy of superconducting gap 
energies for the two bands, with the 7r band becoming 
stronger than the a band at high doping, opposite to the 
order for pure MgB2. In this region, the predicted T c is 
also much lower than observed experimentally, perhaps 
due to uncertainty of the true amount of Al actually in- 
corporated in the superconductor or the inhomogeneous 
distribution of Al. 

In this Letter we present experimental data that vali- 
dates these theoretical predictions better than in earlier 
work, including the inversion of the relative strength of 
superconductivity on the 2 bands when the Al content is 
high. We combine spectroscopy, which tracks the loss of 
hole states on the a band, with heat capacity measure- 
ments, which provides information about the behavior of 
both gaps. We also analyze data in terms of the relative 
change of the unit cell volume v(x) = V(x)/Vo, where 
Vo = V(0) is the unit-cell volume for pure MgB 2 , and 
V(x) is the unit-cell volume obtained from x-ray diffrac- 
tion data after reaction. Since the crystal lattice changes 
systematically with added Al, this procedure gives a bet- 
ter volumetric representation of the Al content actually 
present, on average, in the superconducting phase. Thus, 
we can assert that the behavior that is observed is intrin- 
sic to Al-doped MgB 2 . 

Two sets of samples were prepared using reactions at 
opposite limits of what is routinely used in the literature 
in terms of time and temperature: a short reaction "A" , 
850 °C for 1 hr, which was stopped just short of full con- 
sumption of reactants (thus the superconducting phase 
is "as formed"); and a very hot, long reaction "B", 1200 
°C ramped down to 700 °C over > 80 hours like that in 
|2l|. Full details of the structural characterization will 
be described in more detail elsewhere. It is important 
to note here that the crystal lattice parameters for sam- 
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FIG. 1: Critical temperature as a function of the normalized 
volume of the unit cell taken from x-ray diffraction data for 
samples made by a short reaction (A) and a long, hot reaction 
(B). The curves are guides to the eye. The axis along the top 
edge indicates Al concentration calculated for a simple rule of 
mixtures between MgB2 and AIB2 (Vegards law). 



pie set B generally obeyed a simple mixing rule between 
MgB2 and AIB2 endpoints (Vegards law) and the compo- 
sition measured after reaction by energy-dispersive x-ray 
spectroscopy in a scanning electron microscope. These 
samples thus appeared to be homogeneous in composi- 
tion. Samples from set A, on the other hand, generally 
incorporated Al more slowly than Mg, and also displayed 
other signs of inhomogeneity. The x-ray diffraction char- 
acterization did not indicate the presence of any Mg-Al 
ordering [z| for sample set B, while broadening of the 
(002) peaks was evident for set A for 0.15 < x < 0.25. 

Despite the evidence for inhomogeneity in one sample 
set, superconducting properties analyzed in terms of v for 
both sets of samples fell onto a single curve, as for T c (v) 
shown in fig. ^ Uncertainty bars on this plot represent 
the 10 to 90% width of the T c transition in a 1 mT back- 
ground field, measured after cooling in zero field. The 
T c (v) curve exhibits a steep fall for v > 0.96, down to a 
critical temperature of about 7 K. For v < 0.96 the in- 
ductive transitions clearly showed bulk supercondcutiv- 
ity, where in this region T c extrapolates to at v ss 0.93. 

The large grain size produced by reaction B provided 
excellent samples to probe the electronic structure with 
electron energy loss spectroscopy (EELS) [23. l24|. We 
used a field-emission scanning transmission electron mi- 
croscope (STEM) system used previously to probe the 
electronic structure of individual grains of pure MgB2 
PH . EELS reduces the effects of surface contamination 
because the probe beam passes through the crystal vol- 
ume. Furthermore, special crystal orientations can be 
selected easily in the large thin area of a polycrystalline 
sample. Good correspondence between EELS and XAS 
was obtained in a study of the hole states in pure MgB 2 
|26j . In the present experiment, we extend those analyses 
to Al doped samples from set B. 
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FIG. 2: STEM- EELS data are shown in plot (a). The dashed 
line indicates a Gaussian fit to peaks at higher energy, which 
was subtracted to estimate the pre-peak intensity. The inset 
shows the normalized integrated intensity of the pre-peak at 
186-189 eV. Plot (b) shows simulated EELS spectra obtained 
by density-functional theory using the virtual crystal approx- 
imation. The inset shows the calculated band structure for 
the important a band. In both plots (a) and (b), the Fermi 
energy is fixed and the curves are offset in intensity for clarity. 



The different natures of the electronic p xy states and 
p z states, which are associated with the a band and the 
7T band respectively, lends a key advantage to the spec- 
troscopy experiment. In pure MgB2, the p xy density of 
states (DOS) has a steep edge as a function of energy, 
extending from several eV below Ep to about 0.9 eV 
above Ep pfij . The unoccupied p xy states farther above 
Ep are few in number for about 5 eV, where a sharp 
rise in the p xy DOS marks the main boron K edge. In 
contrast to these sharp features of the a band, the p z 
DOS is essentially constant with energy over that same 
range, providing a relatively constant background to the 
spectroscopy measurements. These differences are am- 
plified further by using a large collection aperture and 
aligning the electron beam normal to the boron planes, 
which produces measured intensity that is dominated by 
the momentum transfer into 
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The essential features of the EELS data, shown in 
fig-Et; are a pre-peak at 186-189 eV energy loss followed 
by a strong series of peaks at 190 eV and above. Here, 
the spectra are aligned at 187 eV, which represents Ep. 
The pre-peak is derived from core electrons being excited 
into the unoccupied p xy states at the top of the a band. 
As Al is added, a systematic loss of pre-peak intensity oc- 
curs, indicating the rise of Ep into the flat portion of the 
a band. This behavior is very similar to that measured 
by XAS 0, which included the intensity due to the p z 
states. The main K edge also shifts downward, which will 
be described in a separate paper. A quantitative assess- 
ment of the p xy DOS in this region as a function of doping 
was made by integrating the pre-peak intensity and nor- 
malizing to the integrated intensity for x — 0. The inset 
of fig. EK shows the 2 different approaches used: First, a 
2 eV wide Gaussian curve was fit to the peak intensity 
at approximately 187.5 eV, the width being due to the 
broadening of the STEM. Second, the tails of main edge 
peaks at ~ 192 and ~ 195 eV were fit by a 6 eV wide 
Gaussian (dashed curves in fig. EK) and subtracted, and 
the remaining intensity was integrated from 185 to 191 
eV. As shown in the inset to fig.Et: extrapolating either 
integrated intensity gives the Al content range where the 
a band DOS becomes too small to contribute appreciable 
intensity, which evidently occurs for x > 0.32. 

The STEM-EELS data are in good agreement with 
density-functional calculations based on the virtual crys- 
tal approximation, shown in fig. Eb- Here, the DOS cal- 
culations have been smoothed to simulate the broadening 
in the experiment. Fig. Eb shows that there is no pre- 
peak intensity for x — 0.30, at which point Ep intersects 
the a band at the T point (fig. Et> inset). The remaining 
flat portion below 191 eV is due to the small amount of 
p z states that are sampled by the simulated collection 
conditions, as indicated by the lack of any change in this 
region between x = 0.30 and x — 0.50. The overall cal- 
culations gave results very similar to those reported in 
|l8|. including the position of Ep relative to the bands 
as shown in the inset of fig. Et>- 

Heat capacity measurements were coordinated with 
the EELS analyses to monitor the superconducting gaps 
for the a and ir bands, A a and A,r, and their DOS. Heat 
capacity data were collected using a commercial system 
and analyzed with a 3-parameter fit (the so-called 'alpha' 
model ^27]) to extract the behavior of A CT and A n as a 
function of v. Generally our data were very similar to 
those reported elsewhere [l7ll28j . 

The results of the alpha-model fitting are shown in 
fig. EJi. These data show a strong decrease of A CT with 
decreasing v (i.e. with higher Al content). The 7r gap 
appears to hold constant for v > 0.98 and then decreases 
for smaller v. However A w remains well below A^ un- 
til their values become equal at v ss 0.96. This merging 
point corresponds to the unit cell volume where the T c 
data changes curvature rather sharply. As with the T c 
data in fig. 2] little difference is seen in the trends for 
samples made by either reaction A or reaction B, sug- 
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FIG. 3: The superconducting gaps for both sample sets ex- 
tracted from heat capacity data are shown in plot (a). The 
inset shows an example of the data obtained for a sample with 
x ~ 0.2. Curves are guides to the eye. The Sommerfeld co- 
efficients extracted from the alpha modef fitting are shown in 
plot (b). Dotted lines are again guides to the eye. 



gesting that these behavior are intrinsic to the average 
Al content in the grains. The merging point of the gaps 
corresponds to x = 0.33 according to Vegard's law. 

For still smaller v, two gaps could not be distinguished, 
even though magnetization measurements demonstrate 
that bulk superconductivity survives, albeit only below 
7 K. These data extrapolate to zero gap for v w 0.93 
(x w 0.55 from Vegard's law). The alpha model also 
provides information about the Sommerfeld coefficients 
for the bands, which are proportional to the individual 
DOS. Hence, the band responsible for this state can be 
inferred. In fig-Et; it can be seen that a progressive loss 
of the a band DOS occurs with decreasing v (increasing 
x), extrapolating to zero for v = 0.96. By contrast, the 
7T band DOS increases with v. The a and it DOS have 
approximately equal values at d k 0.98. These trends 
suggest that it is the tt band that is the stronger band 
for v < 0.96, indicating a reversal of the 2-band hierarchy 
in that regime. 

We can summarize our results in terms of two impor- 
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tant implications for understanding superconductivity in 
doped magnesium diborides. First, our data are gener- 
ally more consistent with the predictions of theory and 
less consistent with the existing experimental studies on 
Al-doped MgB2. In particular, a much stronger rate of 
decrease of T c and A a with composition (here derived 
from Vegard's law) has been found, indicating a stronger 
correlation between the loss of strong superconductivity 
on the a band and the topological transition of the Fermi 
surface |2^| . T c falls from 39 to 7 K by the time Ep cov- 
ers the <7 band at the T point in reciprocal space (at 
x « 0.33), confirming the loss of Fermi-surface area as 
the primary factor for T c in this region. 

Second, for 0.33 < x < 0.55 the system is essentially a 
7r band superconductor with T c values below 7 K. Super- 
conductivity could be induced by the n band onto the a 
band hole pockets, but their rapidly shrinking area makes 
their influence on the overall behavior small. The com- 
plete loss of superconductivity at x « 0.55, when the a 
band is filled completely, correlates with the loss of the 
soft phonon modes in infrared absorption experiments 
[1,|3(J|. Thus, in this region the proximity of the a band 
to Ep still provides interband electron-phonon coupling, 
even though the a band DOS does not strongly influence 
T c . 



Thus, the motion of the flat a band near the Fermi 
surface explains the key characteristics of superconduc- 
tivity in Al doped MgB 2 . In this sense, superconduc- 
tivity in MgB2 can be thought of as being derived from 
stable AIB2, where the ir bonding keeps the position of 
the 7r band relatively fixed while the flat portion of the 
a band is moved up through Ep by the addition of Mg. 
An interesting question, therefore, is whether T c contin- 
ues to increase as the a band moves upward away from 
Ep for Na or Li doping (if it is possible) . Since the soft 
phonon modes should again stiffen as the a band moves 
higher, pure MgB2 may lie close to the maximum of the 
T c vs. electron doping behavior. Alternatively, it should 
be possible to take other diborides with stable ir bonding, 
such as ZrB2, and alloy them with hole donors to drive 
them in the direction of a lattice instability to produce 
superconductivity, at least on the tt band. 
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